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Abstract. The paper is divided into two parts. In the first part we 
lay down the foundation for defining the joint annihilation-preservation- 
creation decomposition of a finite family of, not necessarily commutative 
random variables, and show that this decomposition is essentially unique. 
In the second part we show that any two, not necessarily commutative, 
random variables X and Y, for which the vector space spanned by their an- 
nihilation, preservation, and creation operators equipped with the bracket 
given by the commutator, forms a Lie algebra, are equivalent, up to an 
invertible linear transformation to two independent Meixner random vari- 
ables with mixed preservation operators. In particular if X and Y com- 
mute, then they are equivalent, up to an invertible linear transformation 
to two independent classic Meixner random variables. To show this we 
start with a small technical condition called "non-degeneracy" . 

1. Introduction 

This paper continues the work from [12] . by considering examples in di- 
mensions higher than one. The multi-dimensional case is much more difficult 
than the one-dimensional one. We will restrict our attention to the case 
d = 2. The multi-dimensional case gives us the opportunity to work also in 
the non-commutative framework, by considering "random variables" that do 
not commute. In section 2 we introduce the creation, preservation, and an- 
nihilation {APC) operators for a finite family of not necessarily commutative 
random variables. In section 3, we make the crucial observation that the vec- 
tor space spanned by the identity and APC operators of most of the Meixner 
random variables, equipped with the commutator as a bracket, forms a Lie 
algebra. We use this observation, to define the notion of "Meixner random 
vector (X, Y) of class A4l" ■ We also give two important examples of such 
random vectors, one in which X and Y commute, and another in which they 
do not commute. In section 4, making a simple and natural assumption, that 
we call "non-degeneracy", we classify all Meixner random vectors (X, Y) of 
class Aii, by showing that through an invertible linear transformation, they 
can be reduced to the two examples from Section 3. 
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2. Background 

Let (H, (•, •)) be a Hilbert space over K. Let Xi, X2, . . . , be d symmetric 
densely defined linear operators on H. We denote by A the algebra generated 
by Xi, X2, • • • , Xd- We assume that there exists an element <f> of H, such that 
<j> belongs to the domain of g, for any g G A. We fix such an element 4> and 
call it the vacuum vector. 

Definition 2.1. H^e caW any element g of A, a random variable. For any g 

in A, we define: 

E[g] := {g<PA) (2-1) 

and call the number E[g] the expectation of the random variable g. Finally, 
we call the pair (A, 4>) a probability space supported by H . 

The above definition follows the classic GNS representation (see [13] and 
|15j). We restrict ourselves to finitely generated algebras, but this is not 
necessary. We always work with unital algebras A, that means we assume 
that I £ A, where I denotes the identity operator of H. It is not hard to 
see that, for any g in A, there exists a polynomial p(x\, X2, • • • , Xd), of d 
non-commutative variables xi, x%, . . . , x^, such that g = p(X\, X2, • • • , Xj). 
We also introduce the following equivalence relation: 

Definition 2.2. Let (A, (j)) and (A' , <j>') be two probability spaces supported 
by two Hilbert spaces H and H' , and let E and E' denote their expectations. 
Let X\, X2, ■ ■ ■ , Xd be operators from A, and X[, X' 2 , . . . , X' d operators from 
A' . We say that the vector random variables {X\, X2, ■ ■ ■ , Xd) and (X' 1} X' 2 , 
. . . , X' d ) are moment equal and denote this fact by (X\, X2, ■ ■ ■ , Xd) = (X[, 
X' 2 , ■ ■ ■ , X' d ), if for any polynomial p{x\, X2, ■ ■ ■ , Xd) of d non-commutative 
variables, we have: 

E[p(X l ,X 2 ,..., X d )} =E'[p(X[,X' 2 ,...,X d )]. (2.2) 

For any non-negative integer n, we define the space F n as being the set of 
all vectors of the form p(X\, X 2 , ■ ■ ■ , Xd)4> G H, where p is a polynomial of 
total degree less than or equal to n. It is clear that F n is a finite-dimensional 
subspace of FL. Being finite dimensional, F n is a closed subspace of H, for all 
n > 0. We have: 

F cFiCF 2 C-- - CH 

We define Go := Fq, and for all n > 1, G n := F n Q that means G n is 

the orthogonal complement of F n _\ into F n . For any n > 0, we call G n the 
homogenous chaos space of order n generated by X±, X2, ■ ■ ■ , Xd- We also 
define the space: 
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and call H the chaos space generated by X\, X2, • • • , It is not hard to 
see that H is the closure of the space A4> '■= {g4> \ g G A} in H. 

It is also easy to see, based on the fact that X±, X2, ■ ■ ■ , Xd are symmetric 
operators, that we have the following lemma: 

Lemma 2.3. For any index i £ {1, 2, d} and any non-negative integer 
n: 

XiG n _L Gk, (2.3) 
for all k 7^ n — 1, n, and n + 1. 

See also [1] for a proof. It follows now, mutatis mutandis as in [I], that for 
any i € {1, 2, ... , d}, there exist three operators a~(i), a°(i), and a + (i), called 
the annihilation, preservation, and creation operators, respectively, such that: 

Xi = a-(i) + a°(i) + a + (i). (2.4) 

In (|2.4p . the domain of Xi, a~(i), a (i), and a + {i) is understood to be A(j>. 
It is important to remember that, for any i E {1, 2, . . . , d} and n > 0, 
a~(i) : G n -> G n _i, a°(z) : G n -> G„, and a + (i) : G n -> G„ + i, where 
Gr_i := {0} is the null space. 

If Y" and Z are two operators, then we define their commutator \Y, Z] as: 

[Y, Z] := YZ - ZY. 

It is also not hard to see that the operators X\, X2, ■ ■ ■ , Xj commute among 
themselves if and only if the following conditions hold, for any i, j 6 {1, 2, 
...,d}: 

[a-(i),a-(j)] = 0, (2.5) 

[a-(i),o°(j)] = [a-(j),a°(i)], (2.6) 

[a"(i),a+(j)] - [a-(j),a + (i)] = [a°(j), a°(i)] . (2.7) 

These conditions are derived, as in [T] and [2], from the fact that, for all n > 0: 

pQ, Xj] : G n — > G n -2 © G n _i © G n © G n+ i © G n+ 2, 

by projecting the equality [_Xj,Xj] = on the spaces G n _2, G n _i, and G n , 
respectively. Once we know that these projections are zero, we can conclude 
that the other two projections (on G n+ \ and G n+ 2) are also zero, from du- 
ality arguments, since a? is a symmetric operator, while a + is the adjoint of a~. 

If X\ , X2 , ■ ■ ■ , Xd are classic random variables defined on the same prob- 
ability space (0, J 7 , P) and having finite moments of any order, then we can 
take H = L 2 (Q, T , P) and <j> = 1, i.e., the constant random variable equal to 
1. We then regard X\, X2, ■ ■ ■ , X^ as multiplication operators on the space 
Al C H, where A is the algebra of the random variables of the form p(X\, X2, 
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. . . , Xd), where p is a polynomial of d variables. It is clear that X{Xj = XjXi, 
for all 1 < i < j < d. 

Conditions (I2.5p . (I2.6p . and (I2.7P separate the Commutative Probability 
from the Non-Commutative one. 

Definition 2.4. Let H be a Hilbert space, (A, </>) a probability space sup- 
ported by H, and {^Q}i<i<d elements of A, that are symmetric operators. Let 
{a~.}i<i<d, { a Xi}i<i<d, and {a x .}i<i<d be three families of linear operators, 
defined on subspaces of H, such that, belongs to the domain of a^o^ • • • a^™, 
for all n > 1, (i x , i 2 , ... , i n ) G {1, 2, . . . , d} n , and (e 1 , e 2 , e n ) £ {-, 
0, +}". We say that these families of operators form a joint annihilation- 
preservation-creation (APC) decomposition of {Xi}\<i<d relative to A, if the 
following conditions hold: 

Xi = a-.+a° x .+a+, (2.8) 

[at.)*\M = a-\Atf>, (2-9) 
a~H n C Hn-i, (2.10) 
a° x .H n C H n , (2.11) 

for all 1 < i < d and n > 0, where H-i := {0}, Hq := and is the 
vector space spanned by all vectors of the form a+^aj.^ • • • a+ (f>, where i\, i 2 , 
. . . , ik G {1, 2, . . . , d}, for all k > 1. 

We call a~. an annihilation operator, a x . a preservation operator, and a x . a 
creation operator, for all 1 < i < d. 

For all 1 < i < d, by (a x .)* \A(f> = a~.\A<p we mean: 

(a+u,v) = (u,a-.v), 

for all u and v in A(p, where (•, •) denotes the inner product of H. Since, for 
all i £ {1, 2, . . . , d}, Xi is a symmetric operator, we conclude from (j2.9fl . that: 

{alX\M = a% t \M- (2-12) 

Let us observe that for any X\, X 2l Xd G A, we can consider the 
unital algebra A' C A generated by X\, X 2 , . . . , Xd- Doing the construction 
described before, by considering, for all n > 0, the space F n of all vectors of 
the form p(X\, X 2 , . . . , Xd)4>, where p is a polynomial of degree at most n, 
then defining G n := F n Q F n _i, and so on, we can construct the annihilation, 
preservation, and creation operators a~(i), a°(i), and a + (i) of Xi, respectively, 
where 1 < % < d. It is now clear that these operators form a joint (APC) 
decomposition of {Xj}i<j<d relative to A'. We call this decomposition the 
minimal joint (APC) decomposition of {Xi\\<i<d- We can prove the following 
lemma about the uniqueness of the (APC) decomposition, which justifies the 
choice of the word "minimal" . 
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Lemma 2.5. Let {Xi}\<i<d be a family of symmetric random variables in a 
non- commutative probability space (A, (ft), and {a^.}x<i<d, {a Xi }i<i<d, o-nd 
{ a ti}i<i<d a joint (APC) decomposition of {Xj}i<j<^ relative to A. Let A' 
be the algebra generated by {^Q}i<i<d, and {a~ {i)}\<i<d, {a°(i)}i<i<d, and 
{a + (i)}i<i<d the minimal joint (APC) decomposition o/{Xj}i<j<^. Then for 
any i E {1, 2, ... , d} and any e G {— , 0, +}, we have: 



a% t \A'4> = <f{i)\A'<ft. (2-13) 
Moreover, if A" denotes the algebra generated by uf =1 {a~ + }, then 



A"<P = A'cft (2.14) 

and {a~.}i<i<d, {a Xi }i<i<d, o.nd {a^}i<i<d is a joint (APC) decomposition 
of {Xi}i<i< n relative to A". 

Proof. As before, for all n > 0, let F n be the space of all polynomials of 
d non-commutative variables of degree at most n (in which the variables x%, 
x 2 , . . . , x n , are replaced by Xi, X 2 , ■ ■ ■ , X n ), G n := F n F n _i, and H n the 
space spanned by all vectors of the form a+^aj. • • ■ a + (ft, where i\, i 2 , ■ ■ ■ , 
i n € {1, 2, d}. 

Claim 1. For all n > 0, H n C F n . 

We prove this claim by induction on n. For n = 0, we have Hq = Fq = M.<ft. 
Let us assume that for some n > 1, we have .fffc C for all A; < ra — 1, and 
prove that _ff n C F n . We need to show that d£. a J. • • • (ft £ F n , for all ij., 
i2, i n £ {1, 2, Indeed, we have: 



X^Xi 2 . . . X in c 



(at. + a„. +a„, )(at. -\-a„. +a„. ) • • • (alt. +a„. + a„. )0 

\ ^2-^ ^i-^ / \ i>2 *2 *2 *ra *ra *n 

cdr. </> + a!. 1 at 2 ■ ■ ■ at d <ft. 

■"•ill ' / j ^i-^ %c 2 I'll 



CLrp . CI™ . 



(ei,e2,...,e„)e{-A+} n \{ (+,+,-,+)} 



Because at least one operator from the terms of the last sum is a preservation 
or annihilation operator, it follows from (|2.11[) and (|2.10|) . that: 

(ei,e 2 ,...,en)e{-,0,+}»»\{(+,+,...,+)} 
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Thus, using the induction hypothesis, we get: 



= x h x i2 . . . x in cp - a x h a Z 2 ■ ■ ■ a x d ln 4 

(ei,e 2 ,...,<E„)6{-,0, + } n \{(+, + ,-, + )} 

e F n - (tf + #i + • • • + 
C F n — F n _i 
C F n . 

Claim 2. For all m ^ n, -ff m _L _ff n (as subspaces of H). 

We prove by induction on n, that for all n > 1, we have: -ff n _L if^, for all 
< k < n- 1. 

For n = 1, we need to prove that -f/i _L Hq, which reduces to showing that 
for all 1 < % < d, a^.(j) 1 </). If (•, •) denotes the inner product of H, then it 
follows from condition (|2.9[) and (|2.10p that: 

= (0,0) 
= 0. 

Thus _L (f). 

Let us suppose that H n _L H^, for all < k < n — 1, and prove that 
H n+ \ _L i7 r , for all < r < n. To do this, we need to show that if r < n, then 

for all ii, 12, ■ ■ ■ , i n +i) Ji, J2, • • • , jr € {1, 2, . . . , d}. Using again the condition 
(12. 9p . because i/ ra +i and i/ r are both contained in F n+ \, and F n+ i C .40, we 
get: 

(at. al". • • • at. 4>, at. at. ■ ■ ■ at. (t>) 
= (at. at. ■ ■ ■ at. d>, a~. at. at. ■ ■ ■ at. (£>)■ 

\ *3 n -\- 1 ' 1 Jl J 2 Jr 

From (|2,10p . we know that a~ a+ a£. ■ ■ ■ at. 4> G -Hr-l- Since r < n, we have 
r — 1 < n — 1, and thus, it follows from the induction hypothesis that: 

(ajr. alt. • • • alt. 0, a~ alt. at. • • • at. <t>) = 0. 

\ j,j 2 j,j 3 J '»tt+1 r x Jl ■ L 32 IT ' 

Claim 3. For all n > 0, F n = H Q © H i © • • • © where "ffi" denotes the 
orthogonal sum. 

We know from the previous two claims that Hq © Hi © • • • © H n C F n . On 
the other hand for any monomial X^X^ . . . Xi m (j), of degree m, where m < n, 
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we have: 

x h x i2 . . . x im( j) = J2 a \ a % 2 •••<!„ # 

(ei,e 2 ,...,e m )G{-,0, + } m 

g #o © Hi © • • • e H m . 

Since, these monomials are spanning F n , we get: F n C Hq © Hi © • • • © i? n . 

It follows now from Claim 3, that for all n > 0, 

# n = [iT © Hi © • • • © H n ] G [H © iJi © • • • © # n _i] 
= F n ©F„_! 
= G rt . 

Let n > be fixed and let G H n , then for any 1 < i < d, we have: 

XiW = at.w + a°L.w + a~.w. 

j^i a^i 

Since a+.w G H n+ i = G n+ i, a Xi w G iT n = G n (from ()2.10p ). and a~.w G 
H n _i = G n -i (from (|2.11j) ). we conclude that a x .w = a + (i)w, a x .w = a°(i)w, 
and a~ u; = a~(i)w, and the proof of (12, 13ft is complete. 

Moreover, 

A" $ = U n > (ffifc = o^fc) 
= U n >o-F n 
= A'<f>. 

Let us observe now that, for each 1 < i < d, X; = a~ + a® + a+ G A". 
Thus {^i}i<j<d can be regarded as random variables in the non-commutative 
probability space (A", (f). It is now clear that {a~}i<j<d, {a°.}i<i<d, and 
i a xi}i<i<d is a joint (APC) decomposition of {Xj}i<j< n relative to A". □ 

Using the above notations, we close this section by defining the number 
operator. 

Definition 2.6. We define the number operator N as the linear operator 
whose domain is A'cfi, defined by the formula: 

Nu = nu, (2.15) 

for all n > and all u G H n . 

3. MEIXNER RANDOM VECTORS OF CLASS M L 

In the one-dimensional Meixner random variable X, with infinite 

support, has the Szego-Jacobi parameters: a n = an + ao and u n = j3n 2 + (t — 
f3)n, for all n > 1, where a, f3, and t are real numbers, such that /3 > and 
t > 0. The Meixner random variables (r.v.), with infinite support, are divided, 
up to a re-scaling, into five sub-classes: Gaussian, Poisson, negative binomial, 
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gamma, and two parameter hyperbolic secant r.v.. A Meixner random variable 
X, with finite support, has the following Szego-Jacobi parameters: 



and 



LJ r . 



an + ao if n < k — 1 
if n > k 



(3n 2 + (t- 0)n if n<k — l 
if n > k 



where A: is a natural number equal to the number of different values that X 
takes on with positive probabilities, and a, /3, and t are real numbers, such 
that if k > 2, then t > and t + f3(k - 1) > 0. See @], [9], and [H]. 

Thus we can write the class A4 of all Meixner random variables as: 

M := M u UM b , (3.1) 

where Ai u represents the class of all Meixner random variables with infinite 
(unbounded) support, and A4{, the class of all Meixner random variables with 
finite (compact) support. 

Let us consider now a Meixner random variable with infinite support, having 
the Szego-Jacobi parameters: a n = an + ao and cj„ = (3n 2 + (t — (3)n, for all 
n > 1. This is equivalent to the fact that the commutators of the terms of its 
minimal (APC) decomposition X = a~ + a + a + are: [a~ , a + ] = 2(3J\f + tl, 
and [a~, a ] = era", where A/" denotes the number operator, i.e., the linear 
operator that maps the n— th orthogonal polynomial f n (X), generated by X, 
to nf n (X), for all n > 0. It also follows that a = aj\f + a^I. Thus, if a ^ 
or a = (3 = 0, then we have [a - , a + ] = 2(/3/a)a° + [t — 2ao(/3/a)K and [a~, 
a ] = cra~, where /3/a := if a = /3 = 0. By duality we also obtain; 

[a°,a+] = [(a°T,(a-y} 

,0i* 



= [a-,a L 
= (cra~)* 
= cra + . 

This means that the vector space spanned by the operators a - , a , a + , and I is 
closed with respect to taking the commutator. Thus this vector space equipped 
with the commutator [•, •] forms a Lie algebra. If we define (3 /a := oo, for 
a = and f3 ^ 0, then we can split the class Ai u into two subclasses: M u ,f 
and Ai u ,oo- -M- U j is the class of all Meixner random variables, with infinite 
support, for which a^0ora=/? = 0. -Mu,oo is the class of all Meixner 
random variables, with infinite support, for which a = and (3 ^ 0. Thus an 
element of Ai u ,oo is a random variable having a symmetric (symmetry about 
ao) re-scaled two parameter hyperbolic secant distribution. Moreover, the 
following proposition holds: 
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Proposition 3.1. The Meixner random variables of class M u ,f> are exactly 
those random variables X = a~ +a° + a + , having finite moments of all orders 
and infinite support, for which the vector space W spanned by a~ , a , a + , and 
I, equipped with the commutator [■, ■], forms a Lie algebra, where I denotes 
the identity operator. The Gaussian random variables are the only symmet- 
ric (symmetric about a number) random variables, for which the vector space 
spanned by a" , a + , and I, equipped with the commutator [■, ■], forms a Lie 
algebra. 

Proof. (=>) If X is a Meixner random variable of class M u ,f, then we 
have already explained why the vector space (W, [-, •]) forms a Lie algebra. 

(<=) Let us suppose that (W, [•, •]) is a Lie algebra. Let {f n }n>o be the 
sequence of orthogonal polynomials generated by X, and let G n := E/ n , for 
all n > 0. Since (W, [•,•]) is a Lie algebra, there exist four constants a, c, d, 
and e, such that: 

[a~,a°] = aa~ + ca° + da + + el. (3.2) 

Thus for all n > 0, 

[a-,a°]f n (X) = aa-f n (X) + ca°f n (X) + da + f n (X) + eIf n (X). 

This is equivalent to: 

(otn — a n -i)oj n fn-i(X) = au n f n -i(X) + (ca n + e)f n (X) + df n +i(X), 

for all n > 0, where a_i := 0. Since X has an infinite support, f n -±(X), 
f n (X), and f n+ i(X) are linearly independent, for all n > 1. Because oj n > 0, 
if we look at the coefficient of f n -\ hi the last equality, we get a n — a n _i = a, 
for all n > 1. Thus: 

n 

a n = «o + ^( a k - ctk-i) 
k=l 

n 

= «o + 

k=l 

= «o + an > 

for all n > 1. 

Because (W, [•, •]) is a Lie algebra, there exist four constants p, q, r, and s, 
such that: 

[a~,a + ] = pa~ + qa° + ra + + sL. (3-3) 

Since, for all n > 0, [a~, a + ] and go + sL map G n into G n , pa~ maps G n into 
and ra + maps G n into G n +\, by restricting f|3.3|) to G n , we get: 

[a-,a+]|G n = (qa° + sL)\G n , (3.4) 
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for all n > 0. If we define ojq := 0, then for all n > 0, 

[a",a + ]|G n = (w n+ i - u) n )I\G n 

and 

(ga° + sJ)|G n = (ga n + s)/|G n 

= (qan + qa + s)I\G n . 

Hence, it follows from (|3.4|) . that — oj n = gan + gao + s, for all n > 0. 
Since loq = 0, we obtain: 

n-1 

w n = y^fc+i — 

fc=0 
n-1 

= ^(ga/c + ga + s) 

k=0 

n(n — 1) 
= h (ga + sjn 

= Y n + \ qa ° + s ~ ~2~) n ' 
for all n > 1. Thus X is a Meixner random variable. Moreover, the numbers 
a and /3, of X, satisfy the relation /3 = (q/2)a, which means that if a = 0, 
then /3 = 0, too. Hence X is of class M. u ,f- □ 

We also have the following proposition: 

Proposition 3.2. Let X be a centered random variable, taking on only k dif- 
ferent values with positive probabilities. The space W , spanned by the identity 
and (APC) operators a~ , and a+ of X, equipped with the bracket given 
by the commutator, is a Lie algebra if and only if X is a Meixner random 
variable, taking on only k different values, such that, the constants t and j3, 
involved in its Szegd-Jacobi parameters, satisfy the following condition: 

t = -0(k-l), (3.5) 

and a ^ 0. Conditions 113. 5|) and a ^ are satisfied exactly by the non- 
symmetric binomial distributions, with Krawtchouk orthogonal polynomials. 

Proof. We do the same proof as before, except that due to the fact that 
the last non-zero orthogonal polynomial is fk-i, we have u>k = 0. Thus [a~, 
a + ]fk-i = -Wfc-i/fc-i. Since: 

[o-,a+]|G n _i = {2(/3/a)a° + [t-2ao(/3/a)]/}|G„_ 1 , 

we must have [a~ , a + ]fk-i = {2(/3/a)a° + [t — 2ao(/3 / a)]I} fk-i- This means 

-u)k-i = 2fi(k-l)+t. 

Since w fc -i = P{k - l) 2 + (£-£)(&- 1) it follows that (JOJ) must hold. □ 
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We can split the class of Meixner random variables with finite support as: 

M = M b , f UM b>00 , (3.6) 

where Mbj are the non-symmetric binomials, while Mb )0 o are the symmetric 
ones. 

Let us now define the following class: 

Ml := M u jUM b j (3.7) 

and call it the Meixner-Lie class. We can wrap up our discussion now in the 
following lemma: 

Lemma 3.3. The real vector space spanned by the identity operator, and the 
annihilation, preservation, and creation operators of a random variable, having 
finite moments of all orders, equipped with the bracket given by the commuta- 
tor, forms a Lie algebra if and only if that random variable is of Meixner-Lie 
class. 

It has been known, see for example [6] and [7], that various Meixner classes 
give rise to Lie algebras like: su(l, 1), sl(2), and others. 

We can always assume that X is centered, that means Qo = E[X] = 0, 
since otherwise, we can consider the random variable X' = X — E[X], whose 
w-parameters are the same as those of X, while the a-parameters are a' n = 
a„ - «o, for all n > (in the infinite support case). 

Let us take Lemma [3. 31 as a starting point in defining the notion of Meixner 
probability distributions on M? of class M l . 

Definition 3.4. Let /j, be a probability measure on M 2 having finite moments of 
all orders. Let us denote by (x, y) a generic vector in M 2 . Let a~ , a^., and a% be 
the annihilation, preservation, and creation operators generated by the operator 
X of multiplication by x, and a~ , a?, and a+ the annihilation, preservation, 
and creation operators generated by the operator Y of multiplication by y. We 
call fi a Meixner distribution of class Ml, if the real vector space W spanned 
by the operators a~ , a° x , a~ , a°, Oy , and I, equipped with the bracket [■, •] 
given by the commutator, forms a Lie algebra. 

We extend this definition to the non-commutative case in the following way: 

Definition 3.5. Let {A, cp) be a non-commutative probability space and X 
and Y two random variables from A. Let {a~ , o° , a„ }«e{x,?/} be their minimal 
joint (APC ) decomposition. We say that the pair {X, Y) is a Meixner random 
vector of class Ml if the real vector space W spanned by the operators a~ , a®, 
a%, a~ , a,y, a,y , and I, equipped with the bracket [■, •] given by the commutator, 
forms a Lie algebra. 
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Definition 3.6. Let X and Y be two random variables in a non-commutative 
probability space (A, (ft) supported by the Hilbert space H. We say that the 
random vector (X, Y) is non-degenerate if the vectors X<p, Y<p, and <f> are 
linearly independent in H . In particular if both X and Y are centered, then 
the random vector (X, Y) is non-degenerate if X(j> and Yep are linearly inde- 
pendent. 

Proposition 3.7. If (X , Y) is a non- degenerate random vector, then the 
annihilation operators a~ and a~ , of X and Y , are linearly independent. 

Proof. Let c and d be two real numbers such that 

ca~ + da~ = 0. 
Taking the adjoint in both sides of this equality we get: 

ca x + da^ = 0. 
Since = E[X](p and a? y (\> = E[Y]<f>, we have: 
cX(f> + dYc/) - {cE[X} + dE[Y]) <\> 
= c(a x + a x + a+)(/> + d(a~ + a° + a^)0 — {coP x + da°) <fi 
= 0. 

Thus, because X<j>, Y(p, and eft are linearly independent, we conclude that 
c = d = 0. Therefore, a~ and a~ are linearly independent. □ 

We present now two fundamental examples. 

Example 1. Let X and Y be two independent centered Meixner random 
variables of class Ml defined on the same probability space (0, J 7 , P). Let \i 
be the joint probability distribution of X and Y . We can identify now X and 

Y with the multiplication operators on L 2 (IR 2 , fj,), by the coordinates x and y 
of the generic vector (x, y), respectively. In this way X := a~ + a a x + a x and 

Y := a~ + + a+. Since X and Y are independent, we know from [I], that 
[a^ 1 , a y 2 ] = 0, for all (ei, 62) € { — , 0, +} 2 . Moreover, one can see that: 

[a~,a+] G M/ + lRa°, (3.8) 

[a~,a+] G M/ + lRa°, (3.9) 

[a-,a° x ] G Ra~, (3.10) 

[a° x ,at] G Ma+, (3.11) 

[a-,a° y ] G Ra~, (3.12) 

[a°,a+] G Ma+. (3.13) 

Hence (W, [•,•]) is a Lie algebra, where W is the real vector space spanned 
by /, a~, a x , a x , a~ , a®, and a+. Thus pi is a Meixner distribution on M 2 of 
class Ml, or equivalently (X, 1") is a commutative Meixner random vector of 
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class M.l- 

Example 2. Let T and Z be two independent centered Meixner random 
variables of class Ml defined on the same probability space (Q, J-, P), hav- 
ing the same numbers a = 1, t' = 1 (we are using t' instead of t, since the 
letter t we will be used later as a subscript), and the f3 numbers as follows: 
Pt '■= (l/2)(cp + dr) and Pz '■= (l/2)(js' + kv), where c, p, d, r, j, s', k, and 
v are some given real numbers (the choice of these letters will become more 
clear in the next section), such that: 

cs + dv = 

and 

jp + kr = 0. 

In other words, using the usual dot product "•" on M 2 , and the orthogonal 
relation "_L" generated by it, we have: 

Pt = ~(c,d)-(p,r), (3.14) 

Pz = l(j,k) (3.15) 

(c,d) _L (s',v), (3.16) 
(j,k) J- (p,r). (3.17) 

Relations (|3.14p - (l3,17p are important, for the following reason. Pt and Pz 
must be non-negative. We can make sure that they are not negative by re- 
membering that, the dot product of two vectors is the product of the length 
of the vectors times the cosine of the angle formed by the vectors. Thus we 
can start with two non-zero vectors (c, d) and (p, r) forming an acute angle. 
Then we rotate the semi-lines supporting these vectors by 90°, and choose 
any two non-zero vectors (s' , v) and (j, k) on the rotated semi-lines. In this 
way we know for sure that Pt and Pz are strictly positive and their sides are 
perpendicular, as required by conditions (|3.16p and (|3.17p . 

Let [i be the joint probability distribution of T and Z. As before, we can 
identify now T and Z with the multiplication operators (densely defined) on 
H = L 2 (R 2 , fi), by the coordinates t and z of the generic vector (t, z). In this 
way T = + a® + af and Z = a~ + a® + af[ . Let us consider the following 
symmetric operators: 



and 



X := a t + (pa® + s'a° z ) + 



Y := a~ + (ra° + va° z ) + a+ . 



Let A be the unital algebra generated by a t , a° , af , a z , a° z} and af . Let 
4> := 1 (the constant random variable, of H, equal to 1). Then (A, 1) is a 
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non-commutative probability space supported by H, and X and Y are ran- 
dom variables in (A, 1). 

Claim 1. (X, Y) is a Meixner vector of class Ml- 

Indeed, a joint {APC) decomposition of X and Y, relative to A is: 

a~ = at, (3.18) 

a° x = pa° t +s'a° z , (3.19) 

4 = 4, (3-20) 

dy = a~, (3.21) 

a° y = ra° t +va° z , (3.22) 
a+ = a+. (3.23) 

Using now Lemma [2 .51 we conclude that this (APC) decomposition, restricted 
to the space A'l, where A' is the unital algebra generated by X and Y, is the 
minimal (APC) decomposition. 



Since T and Z are linearly independent we have: 



and 



[a° x ,a° y ] = 0. 



Because T is a centered Meixner random variable of class Ml, we have: 

[a~,a+] = [at, af] 

= ^a° t +t'l 
a 

= (cp + dr)a^ + /. 

Multiplying both sides of the relation a x = pa® + s'a° z by c and both sides 
of a® = ra^ + va a z by d, and then adding the two relations together, since 
cs' + dv = 0, we get: 

(cp + dr)a® = ca x + da®. 

Thus, we get: 

[o-,a+] = ca° x + da° y + I £ W. 

Similarly, we have: 

[a~,a+] = ja° x + ka° y + I £ W. 
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We also have: 

[at,pa° t +s'a° z ] 

P [<h><$] + s ' [ a 7, a °z] 
paa^ + s'(0) 

= P a t 

= pa~ G W. 

Similarly, we have: 

[a-,a° y ] = ra~ £ W, 

[ay,a° x ] = s'a'eW, 

[o-,oJ] = va'eW. 

By duality (taking the adjoint) we also get that the commutators between the 
preservation and creation operators belong to W. 

Thus the space W, spanned by /, a~, a°, a+, a~ , a°, and a+, equipped with 
the commutator [•, •] is a Lie algebra. Hence, (X, Y) is a Meixner random 
vector of class M. l ■ 

Claim 2. If T and Z are both not almost surely equal to 0, then the centered 
random vector (X, Y) is non-degenerate. 

Indeed, in this case the vacuum vector is <p = 1. Because a°l = a^l = 
(since T and Z were assumed to be centered from the beginning), we have 
XI = Tl = T and Y\ = Zl = Z. Since T and Z are independent, they must 
be linearly independent, because if T = XZ, for some A / 0, then: 

= E[T]E[Z] 
= E[TZ] 
= E[XZ 2 } 
= XE[Z 2 ] 
+ 0. 

We get a contradiction. Thus XI and Y\ are linearly independent and so, the 
centered random vector (X, Y) is non-degenerate. 

Let us observe that X and Y commute if and only if [a~ , a°] = [a~, a°] ( (I2.6P 
is the only axiom of the commutative probability that is not guaranteed to 
hold in this example), which means r = s' = 0, and hence: X = +pa^ + af 
and Y = a~ + va® + a[ are two independent centered Meixner random vari- 
ables, with parameters ax = P, ay = v, fix = Pt, Py = Pzi and t' = 1. 
Moreover, X is of class A^l, since if ax = 0, then f3x = (l/2)(cp + dr) = 0, 
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because p = a x = and r = 0. Similarly Y is of class Ml- Therefore, if 
XY = YX, then Example 2 reduces to Example 1. 

Since the annihilation and creation operators of X and Y are the same as the 
annihilation and creation operators of T and Z, and the preservation operators 
of X and Y are superpositions of the preservation operators of T and Z, and 
because T and Z are independent, we call X and Y, from Example 2, two 

independent Meixner random variables of class M l with mixed preservation 
operators. 

Observation 3.8. If {X , Y) is a non-degenerate Meixner random vector of 
class Ml, then for any invertible linear transformation T : M 2 — > M 2 , (X' , 
Y') := T(X, Y) is also a non-degenerate Meixner random vector of class Ml- 
In particular, if XY = YX, then X'Y' = Y'X'. 

4. The classification of non-degenerate centered 
Meixner random vectors of class Ml 

In this section we show that every non-degenerate centered Meixner ran- 
dom vector (X, Y), of class Ml, can be reduced, through an invertible linear 
transformation, to a random vector equivalent to a vector of two indepen- 
dent centered Meixner random variables X' and Y' , of class Ml, with mixed 
preservation operators. In particular if X and Y commute, then X' and Y' 
are independent. 

We present first the following lemma. 

Lemma 4.1. // {X , Y) is a Meixner random vector of class Ml, then the 
following relations hold: 

(1) There exist some real numbers b, c, d, e, f, g, h, j, k, p, q, r, s, r' , 
s' , u, and v, such that: 



\a x 


,at] 


= bl + ca x + day, 


(4.24) 


a x 


,a+] 


= el + fa x + ga° y , 


(4.25) 


a~. 


,a+] 


= el + fa x + ga° y , 


(4.26) 


.V 


,a+] 


= hi + ja° x + ka° y , 


(4.27) 




' a x] 


= Pa x +qa~, 


(4.28) 




4] 


= Pat+qa^, 


(4.29) 




.<$] 


= ra~+say, 


(4.30) 




4] 


= rat + say, 


(4.31) 
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ciy , a x ] 


= r a x + s'ciy , 


(4.32) 


' +1 
a xi a y J 


= r d x + s dy , 


(4.33) 


a- , aj] 


= ua~+va~, 


(4.34) 


a y , a,y ] 


= ud x + vciy. 


(4.35) 



(2) 



a x i a y ] 


= o, 


(4.36) 


a x ) a y ] 


= o, 


(4.37) 




= 0. 


(4.38) 



(3) If XY = YX, then we can take r' = r and s' = s. 
Proof. Let {G n } n >o be the homogenous chaos spaces generated by X and 

Y. 

1. Since, for each n > 0, the operators [a~,a+], [a~,a+], and [a~,a+] arc 
mapping G n into G n , and because (W, [•, •]) is a Lie algebra, each of these 
commutators must be a linear combination of only three operators from the 
set {I, d x , d x , d x , dy , dy, dy}, namely, /, d x , and dy. Similarly since, for 
each n > 0, [a~,a°], [a~,a°], [a~,a x ], and [a~,a°] map G n into G n -i> each 
of these commutators must be a linear combination of d~ and a~. Therefore, 
there must exist some real numbers 6, c, d, e, /, /i, j, fc, p, q, r, s, r', s', u, 
and w such that: 

[a~,a x ] = bl + cd x + ddy, 
[a~,a+] = el + fd x + gd° y , 
[ a y> a y] = hl + ja x + kay, 
[ a x, a °x] = Pax+ aa y, 

[dy,d° x ] = r'd x +s'dy, 
[ a y' a y] = ud~+vdy. 



Taking the adjoint in both sides of (I4T25D . dOgj) . (OOP . (Q21) . and (|4T34"1) we 
get dH26]), IfOSjL (OP . IfOSjL and respectively. 

2. Since, for all n > 0, [o+, a+] maps G n into G n+ 2, and no operator from the 
set {I, a~, a x , a+, a~, a°, a+}, spanning W, has this property, we have: 

[a+,a+] = 0. 
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Similarly, we can see that: 

[a~,a~] = 0. 

Since, for all n > 0, [a x , a°] maps G n into G n , there are three real numbers a, 
(3, and 7 such that: 

[a x ,a° y ] = al + /3a° x + 70°. (4.39) 
The left-hand side of (|4.39p is an antisymmetric operator, due to the fact that: 

[a x , a y ] \ a yi a x\ 

= -[a° x ,a% 

while the right-hand side is a symmetric operator. Thus: 

[a° x ,a° y ) = 0. 

3. If XY = YX, it follows from (|2.6[) that [a~ , a°] = [a~, a°], and thus we 
can take r' = r and s' = s. □ 



Let us not forget that the random variables X and Y are assumed to be 
symmetric operators. We are presenting now two important lemmas. 

Lemma 4.2. If (X , Y) is a centered non-degenerate Meixner random vector 
of class A4l, then there exists an invertible linear transformation T% : R 2 — > 
M 2 , such that the coefficients b' , e' , and h! of the Meixner random vector of 
class Ml, (X' , Y') := T(X, Y), from Lemma are b' = h! = 1 and 
\e'\ < 1. 

Proof. Since both X and Y are centered we have a x 4> = E[X]<p = and 
aPy(j) = E[Y]4> = 0. Thus, we get: 

E[X 2 } = {{at + al + a-)X<PA) 

= (X<f>, a-<t>) + (X<f>, a° x <P) + (a-X<f>, 0) 
= O + O + (a-(a-+a°+a+)0,0) 

= {a>ia-x<t>,<l>) + (i a x , a t}4>,0) 
= + ({bl + ca° x + da° y )(f>,(f>) 
= b, 

where (•, •) denotes the inner product of the Hilbert space H supporting the 
non-commutative probability space (^4, eft), in which X and Y are random 
variables. Similarly, we have: 

E[Y 2 ] = h 

and 

E[XY] = <[o-a+]l,l> 
= e. 
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Since (X, Y) is non-degenerate we have X<fi ^ and Yep ^ 0. Because X is 
a symmetric operator, we have: 

b = E[X 2 } 

= {x 2 4>A) 

= {X^X(j>) 

= II X( p II 2 
> o, 

where || • || denotes the norm of H. Similarly, we have h > 0. Let us consider 
the invertible linear change of variables (in fact a simple re-scaling) : 

It is easy to see now that the coefficients b' and h! of (X' , Y') are b 1 = h! = 1. 
Thus by this re-scaling we may assume that b = h = 1. It follows now from 
the Schwarz' inequality, that: 

\e\ = \E[XY}\ 

= KXY&W 

= \(Y<t>,X<t>}\ 

< || Xcf> || • || Y<j) || 

= Vb-Vh 

= 1. 

Thus |e| < 1. We cannot have |e| = 1, since if this were true, then we would 
have equality is the Schwarz' inequality that we used, which would imply that 
Xcj) and Yep are linearly dependent, contradicting the fact that (X, Y) is 
non-degenerate. Hence |e| < 1. □ 

Lemma 4.3. Let (X , Y) be a non-degenerate centered Meixner random vector 
of class Ml- There exists an invertible 2x2 matrix T, with real entries such 
that for the non- degenerate centered Meixner random vector (X\, Y\) := T(X , 
Y) of class M.l, the coefficients q, s, r' , and u, from Lemma\J^j\ o,re q = s = 
r' = u = 0. Moreover, b = h = 1 and \e\ < 1. 

Proof. From Lemma [4.2l we may assume that: b = h = 1 and |e| < 1. Let's 
find first some relations between the values of the coefficients from Lemma [4. 11 
that hold for all non-degenerate Meixner random vectors (X, Y). 

From the Jacobi identity: 

[a-x, [a° x ,ay]} + [a° x , [a° y ,a~]} + [aj, [a~ ,a° x }} = 0, 
using the fact that a~ and a~ are linearly independent, we get: 

sr' = uq (4.40) 
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and 

s(s'-p) = -q(r-v). (4-41) 
From the Jacobi identity: 

[a y ,[a x ,a y ]} + [a x ,[a° y ,a y }} + [a y ,[a-,a x }} = 0, 

we obtain: 

u(s'-p) = -r'(r-v). (4.42) 
Let 7 := s' — p and 5 := r — v. Relations (|4.4ip and ()4.42p become now: 

S7 = -q5 (4.43) 

and 

try = -r'8. (4.44) 
Formulas (fOgJ) . (|4T30l) . (|4T32l) . and (|4TMD can now be written as: 



[ax, a °x] = Pax+Q a y 

[a~,a°] = (v + 5)a- + sa y 

[a y ,a° x ] = f 'a~ + (p + j)a~ 

[a y ,a° y ] = ua-+va y . 

Case 1. If q = s = r' = u = 0, we have nothing to prove and are done. 

Case 2. Let us assume that at least one of the numbers q, s, r' , and u 
is not equal to zero. 

We try now to find two random variables Z w = a w X + f3 w Y, w S {1, 2}, 
where a w and j3 w are real numbers, such that there exist some real constants 
X w and fj, w for which the following relations hold: 

[a~ Zw ,a x ] = Aa~, (4.45) 

[ a i' a °] = / xa ^' (4-46) 

and the matrix: 

Ai fii 

A2 fJ>2 

is invertible. We drop for the moment the index w. It is easy to see that if 
Z = aX+fiY, then an (APC) decomposition of Z is given by: a\ = aa e x +(3a y , 
for all e G {-, 0, +}. Thus: 

= a [a~,al] + /3 [a~,a x ] 

= (ap + f3r')a~ + [aq + fi(p + j)]a y . 
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It follows from here that if we want (|4.45p to hold, then the coefficients ap+Pr' 
and aq + P(p + 7) must be proportional to a and j3. Since we also want (a, 
P) 7^ (0, 0), this fact is equivalent to: 

ap + Pr' aq + P(p + j) 

a (3 

This means the following relation must hold: 

r'p 2 - 7 /3a - qa 2 = 0. (4.47) 

A similar calculation shows that if we want relation (|4.46|) to hold, then we 
must have: 

uf3 2 + 6Pa - sa 2 = 0. (4.48) 

Let us look now at the matrix formed by the coefficient vectors (r', —7, —q) 
and (u, (5, — s) of the unknown vector (p 2 , Pa, a 2 ) from equations (14. 47ft and 
(|4.48j) . This matrix is: 

(t> -7 -q\ 
V u — s J 

Amazingly, it follows from relations (|4.40p . (|4.43p . and (|4.44|) . that the deter- 
minant of any 2x2 sub-matrix of this matrix is zero. Thus the two rows of 
this matrix are linearly dependent and so the equations (|4.47p and (|4.48|) are 
equivalent, unless one of them is the trivial equation = 0. 



Since at least one of the numbers q, s, r', and u is not zero, we know for 
sure that at least one of the equations (|4.47j) and (|4.48j) is quadratic in either 
P or a, and the solution(s) of that equation will also be solution(s) of the 
other one. Let us assume r' 7^ 0, and focus on equation (I4.47p . Let us choose 
a := 1. Thus this equation becomes: 

r'p 2 --fP-q = 0. (4.49) 

The solutions of this equation are: 



Pi = W 9 7 , + 4r/g (4.50) 
2r 



and 



P, = Wg±jgg , (4.51) 
2r 

We need P\ and pi to be real and distinct. If we show this, then the matrix: 

(4.52) 



(Oil «2 > 








V Pi Pi) 




( Pi 


P2 ) 



will have a non-zero determinant, and thus the linear transformation (X, 
Y) 1 — y (Zi, Z2) will be invertible. 
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Thus to achieve our goal we need now to show that the discriminant: 

A = 7 2 + Ar'q 

of the quadratic equation (|4.49p is strictly positive. To show this we use again 
one of the Jacobi identities, applied to the vacuum vector <j>, namely: 

[a~, [a°,o£|] 0+ [a°, [a£,a~]] </)+ [a+, [a~,a°]] <j> = 0. 
Since a® (ft = o^<ft = 0j we get: 

r'-q = -erf. (4.53) 

Let us remember, that from the Schwarz' inequality and the fact that (X, Y) 
is non-degenerate we know that |e| < 1. Thus we have: 

-Ar'q < (r'-q) 2 (4.54) 

2 2 

= e 7 

< I7 2 . (4.55) 

It follows now that A = 7 2 + Ar'q > 0. This inequality must be strict, since if 
A = 0, then we must have equality in (|4.54p . which means (r' + q) 2 = 0. Thus 
we would have q = —r', which would imply: ej = —2r' 7^ 0. Hence 7 7^ 0, 
and thus (|4.55p is a strict inequality, which shows that A > 0. 

It follows now easily from (|4.45|> and (14.460 that, for all (i, j) G {1, 2} 2 , we 
have: 



G Ra~. (4.56) 



Re-scaling the random variables Z\ and Z2, we may assume that their 
coefficients b and h are both equal to 1. □ 

Lemma 4.4. // (X, Y) is a centered non-degenerate Meixner random vector 
of class Ml, whose constants from Lemma \4-1\ satisfy the conditions b = h = 1 
and q = s = r' = u = 0, then: 

cs' + dv = 0, (4.57) 
jp + kr = 0, (4.58) 
fp + gr = 0, (4.59) 
fs' + gv = 0. (4.60) 

Proof. From the Jacobi identity: 

[a~, [a~,o+]] + [o~, [a+,a~]] + [a+, [a~,o~]] = 

we get: 

= [a~, [ay, a^)] + [a~, [a£,a~]] 

= [a~ , I + fa° x + a^] — [a y , I + ca° x + da®] 
= (fp + gr)a~ ~{cs' + dv)a~. 
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Since (X, Y) is non-degenerate, a x and linearly independent. Thus 

we obtain that relations (|4.57p and (|4.59p hold. 

Similarly, from the Jacobi identity: 

[ay, [a~,a+]] + [a~, [ay,a~ ]] + [a+, [a~,a~]] = 0, 
we conclude that (|4.58j) and (|4.60|> hold. □. 

We are ready now to present the main theorem. 

Theorem 4.5. If (X , Y) is a non-degenerate centered Meixner random vec- 
tor, then there exists an invertible linear transformation S : M 2 — > M 2 , such 
that the random vector (X' , Y') := S(X, Y) is equivalent (moment equal) 
to a random vector of two independent Meixner random variables with mixed 
preservation operators of class Ml- Equivalently, the vector space spanned by 
the identity operator and the joint (APC ) operators of X and Y is isomorphic, 
as a Lie algebra, to the vector space spanned by the identity operator and the 
joint (APC) operators of two independent Meixner random variables of class 
Ml, with mixed preservation operators. 

Proof. From the previous lemmas we may assume that q = s = r' = u = 
and b = h = 1. 

From the Jacobi identity: 

[ay, [a° x ,at]] + [a° x , [a+ ,a~]] + [a+ , [a~ ,a°]] = 0, 

we get: 

[a~,pa x ] + + [a+, s'a~] = 0. 

This is equivalent to: 

7 [a~ a+] = 0. (4.61) 

Similarly, from 

[a~, [a° y ,a+]] + [aj, [af,a~]] + [a+, [aj.aj]] = 0, 

we obtain: 

S[a~,a+] = 0. (4.62) 

Since [o~,a^] = [a~,a^"], we conclude from (I4.6ip and (14.621) . that if 7 / 
or 8 7^ 0, then: 

[a',a+] = 0, (4.63) 
and thus we can take e = / = g = 0in (|4.25p . We analyze now two cases: 

Case 1. If 7 ^ or S ^ 0, then [a~, a+] = [a~, a+] = 0. We know 
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from the previous lemma that relations (|4.57p and (I4.58P hold. These are ex- 
actly the orthogonality relations from Example 2, from the previous section. 
Moreover the commutators of the joint (APC) operators of X and Y are ex- 
pressed in terms of / and the joint (APC) operators of X and Y in exactly 
the same way as in Example 1. So we are done. 

Case 2. If 7 = S = 0, then p = s' and v = r. It follows now from the 
previous lemma that: 

cp + dv = 0, (4.64) 
jp + dv = 0, (4.65) 
fp + dv = 0. (4.66) 

Since p = s', [a~ , a x ] = pa~ , and [a~ , a x ] = s'a~, we conclude that, for all 
z E {x, y}, 

[a~,a°] = pa~. 
By duality it follows now that: 

[a° x ,a+] = pat, (4.67) 
for all z E {x, y}. Similarly, we have: 

[a° y ,at] = vat, (4.68) 

for all z £ {x, y}. 

Claim. a x = pj\f and = vM, where J\f denotes the number operator. 

Indeed, for any n > and any (zi, z%, ... , z n ) £ {x, y} n , using the product 
rule for commutators and the fact that a x (p = 0, we have: 



2(< 

° x <P+[c 



[a' a' ■ ■ ■ a' 1 a"0 + la?., at at, • • • at 



z n ) "'xY 1 L X ' z l z 2 

a? 



fc=l 

n 

Z)4r • • a i-i (^ a i) a i + i • • • a t<t> 



/c=i 



k=l 



= P™« 2 •••aJ„<P- 

Hence, for all n > and all £ £ H n , we have = pn^. Since we are inter- 
ested only on the action of a x on the space A'4>, where A' is the unital algebra 
generated by X and Y, and A'cj) = U n >o ©fc =0 Hfc, we conclude that a x = pN '. 
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^ - vAT. 
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Similarly, we can see that 



It follows from this claim and the relation (|4.64|) that: 



I + ca x + da u y 



= I + (cp + dv)Af 
= I + OAT 
= I. 

Similarly, it follows from relations (|4.65p and (14. 66ft that: 

[a-,a° y ] = el 

and 

[a-,a° y ] = I. 

Since |e| < 1, we can make the following invertible linear change of variable: 

( I 1 



(X',Y>) = 
It easy to see now that: 



[ a x» a t>] 



a x» a y> 



[a-,,a° yl ] 



a , , a x , 



a y'i a y> 



(X + Y), 

I, 
■ 0, 

p + v 



:{X-Y) 



(4.69) 



p — v 
p + v 
p — V 



The proof is now complete for the following reason. Since the Lie algebra 
generated by the joint (APC) operators of {X' , 1"') is isomorphic to the Lie 
algebra of the (APC) operators of two independent Meixner random variables 
with mixed preservation operators, and because the joint moments can be re- 
covered from the commutators, as shown in [12], in the commutative case, and 
can easily be extended to the non-commutative case, we conclude that (X 1 , 
Y') is moment equal to a random vector whose components are independent 
Meixner random variables of class Ml with mixed preservation operators. □ 
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Due to the fact that in Example 2, the construction starts with two indepen- 
dent Meixner random vectors T and Z, we can make the following observation. 

Observation 4.6. The Lie algebra generated by the joint (APC) operators 
of a two dimensional non-degenerate Meixner random vector (X, Y) of class 
Ml is isomorphic to a subalgebra o/sl(2) © sl(2) © MI. 

Corollary 4.7. Let n be a Meixner probability measure of class Ml on R 2 , 
such that, \i is not supported by any line of equation ax + by = c, with a 2 + b 2 > 
0. Then, the following statements are true: 

(1) Up to an invertible affine transformation, fi is a product of two Meixner 
probability measures of class Ml on R. That means, there exist an 
invertible linear transformation S : R 2 — > R 2 , a vector c = (ct, C2) £ 
R 2 , and two Meixner probability distributions /ii and ^2 of class Ml 
on R, such that the measure: 

u{B) := ii{SB + c), (4.70) 

can be written as: 

v = ^i©^2, (4-71) 

where SB + c := {S(x, y) + c | (x, y) G B}, for all Borel subsets B of 
R 2 . 

(2) If /j, is not supported by any finite union of lines of equation ax+by = c, 
with a 2 + b 2 > 0, then up to an invertible affine transformation fj, is a 
product of two Meixner probability distributions with infinite support, 
of class Ml on R. 

Final Comment We left out of our discussion the symmetric two parameter 
hyperbolic distributions. The Lie Algebra W that we used in this paper, can 
be enlarged, so that we can include also these distributions and characterize 
the entire Meixner class using this new Lie Algebra. We hope to do this in an- 
other paper. We would like to mention that the recursive relation among the 
orthogonal polynomials, in d variables, that appears in leads to commu- 
tation relations among the annihilation, preservation, and creation operators, 
that cannot be related to the present work, but to the next paper. 
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